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Abstract. In this paper we define the continuous Level Lines Shortening evolution of a 
two-dimensional image as the Curve Shortening operator acting simultaneously and inde- 
pendently on all the level lines of the initial data, and show that it computes a viscosity 
solution for the mean curvature motion. This provides an exact analytical framework for 
its numerical implementation, which runs on line on any image at http://www.ipol.im/. 
Analogous results hold for its afiine variant version, the Level Lines Affine Shortening. 
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1. Introduction 

In [S], [9] there were introduced two new image processing algorithms. Level Lines Short- 
ening and its affine variant Level Lines Affine Shortening simulating a contrast invariant and 
scale invariant evolution of an image by mean curvature motion, respectively affine curvature 
motion. The aim of this work is to rigorously justify that the Level Lines Shortening algo- 
rithm proposed in these works computes explicitly a viscosity solution for the Mean Curvature 
Motion 

du 

(1.1) — = \Du\curv{u), 

respectively that the Level Lines Affine Shortening algorithm provides a viscosity solution for 
the Affine Curvature Motion 

(1.2) ^ = \Du\[curv{u)f''^. 

These equations are particularly interesting since they can be axiomatically obtained from 
the image multiscale theories, as the unique partial differential equations satisfying the most 
desired invariance properties in computer vision. This axiomatic characterization was given 
by Alvarez, Guichard, Lions and Morel in [1]. Caselles et al. [1] realized the potential of 
processing directly the image level lines. They proposed to perform a contrast invariant image 
analysis directly on the set of level lines, or topographic map. A fast algorithm computing the 
topographic map was developed by Monasse and Guichard in [21]. 
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The Level Lines Shortening numerical chain stands for the simultaneous and independent 
curvature evolution of all the level lines for a given function. Curve smoothing by the intrinsic 
heat equation, also called Curve Shortening 



was one of the first versions of curve analysis proposed by Mackworth and Mokhtarian in 
[18]. Here k(x) denotes the curvature vector at x, defined as the second derivative x"{s) with 
respect to any length parameter s. By this (nonlinear) evolution a curve instantly becomes 
smooth, shrinks asymptotically to a circle and develops no singularities or self-crossings. 
Rigorous proofs were given by Gage and Hamilton for convex Jordan curves [13] and later 
extended to embedded curves by Grayson [l^. The Affine Shortening equation 



is a surprising variant of curve shortening introduced by Sapiro and Tannenbaum in [23], |24j . 
Angenent, Sapiro and Tannenbaum [25] gave the existence and uniqueness proofs for affine 
shortening and showed a result similar to Grayson's theorem: a shape eventually becomes 
convex and thereafter evolves towards an ellipse before collapsing. A remarkably fast and 
geometric algorithm for affine shortening was given by Moisan in j20j . 

On the other hand, Osher-Sethian defined and studied in [22] the level set method for the 
motion of fronts by (mean) curvature, for which Chen-Giga-Goto [6] and Evans-Spruck [11] 
provided rigorous justifications. In this setting, the initial curve Sq is considered as the zero 
level set of some function uq. The basic result asserts that the zero level set of the evolved 
function 



does not depend on the choice of the initial data and therefore the evolution is purely geomet- 
rical. Their arguments are based on the notion of viscosity solution of Crandall and Lions, 
that allows one to give a suitable meaning to the (MCM) and (ACM) equations, in the class 
of uniformly continuous functions. We refer to the 'user's guide' of Crandall, Ishii and Lions 
[lUj for further details about viscosity solutions. A mathematical link between the median 
filter and the motion by mean curvature was conjectured by Merriman, Bence and Osher [19] 
and later proved by Barles and Georgelin in [2] using viscosity methods and by Evans in [12] 
using a nonlinear semigroup approach. 

Evans and Spruck checked in [11] the consistency of the level set approach with the classical 
motion by mean curvature. More precisely, they showed that the mean curvature motion 
agrees with the classical motion, if and as long as the latter exists. The results applies 
for a smooth hypersurface, given as the connected boundary of a bounded open set. Their 
arguments are based on comparison techniques with lower barriers for the approximated 
mean curvature motion and strongly use the fact that the hypersurface is the zero level 
set of its (signed) distance function. However, the result does not describe the complete 
behavior of all the level lines of a Lipschitz function. Namely, if we are given a Lipschitz 



(1.3) 




(1.4) 




St = {n(-,t) = 0} 
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function uq which evolves by mean curvature in the viscosity sense, are all of its level lines 
evolving independently by curve shortening? For dimension n > 3 the result is not true, since 
hypersurfaces can develop singularities and possible change topology. Thanks to Grayson's 
theorem, the 2D case has a very peculiar structure which we will take advantage of and show 
that evolving independently and simultaneously by curve shortening all the level lines of a 
function is equivalent to applying directly a mean curvature motion to the functions itself. 

The Level Lines Shortening builds on the above mentioned contributions and connects 
explicitly the geometric approach for curve shortening evolutions and the viscosity framework 
for curvature motions. More precisely, this operator first extracts all the level lines of an image, 
then it independently and simultaneously smooths all of its level lines by curve shortening 
(CS) (respectively affine shortening (AS)) and eventually reconstructs, at each step, a new 
image from the evolved level lines. The chains are based on a topological structure, the 
inclusion tree of level lines as a full and non-redundant representation of an image |5], and 
on a topological property, the monotonicity of curve shortening with respect to inclusion. 
Therefore, the hierarchy of the level lines is maintained while performing the smoothing. 

We show in this paper that the image reconstructed from the evolved level lines is a viscosity 
solution of the mean curvature motion (MCM) (resp. affine curvature motion (ACM)). The 
initial image will be considered as an element of a particular space of functions V5(r2) that 
we term space of very simple functions. This corresponds to bilinearly interpolated images 
defined on a rectangle whose topographic maps contain only Jordan curves. The set of 
very simple functions arises naturally in image processing, since level lines corresponding to 
noncritical levels are sufficient to grant an exact reconstruction of the digital image. 

In this way, the described algorithm corresponds exactly to its numerical implementation 
[7] and has the advantage of satisfying both numerically and analytically all the invariance 
properties required by the scale space in question. 

The paper is organized as follows. In section §2 we define the class of very simple functions 
as approximations for Lipschitz functions. Section §3 is devoted to the definition of the Level 
Lines Shortening evolution, as an operator acting both on crowns of Jordan curves and fiat 
areas, and on very simple functions. In the last section §4 we give the equivalence result. The 
result is then extended to general Lipschitz functions, by first approximating the initial data 
by very simple functions and then using standard stability properties of viscosity solutions. 

2. Modeling Level Lines Shortening 

2.1. Crowns of Jordan Curves. A Jordan curve is a one to one continuous map from the 
unit circle into M?. A Jordan curve S splits the plane in two connected components. 
We denote by Lnt{T,) the open bounded component and by Ext(Ti) the open unbounded 
component. 

Definition 2.1 (Partial order). Let and he two Jordan curves. We say that Y} 
surrounds (strictly) and we write T} < Tp' (Y} -< iP) if Int{T}) C Int{Y?') (respectively 
IntiY^) C Int{T,'^)). 
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This defines a partial order on the set of planar Jordan curves. 

Definition 2.2. We say that a Jordan S curve is piecewise C^, or in Cp'^{S^) if it has 
finite length /(S) > 0, any length parametrization is piecewise and if at each discontinuity 
point for the tangent there are left and right tangent vectors, which are not collinear. 

Definition 2.3. We say that a sequence of curves S" converges in Cp''^{S^) to a curve 
S if, denoting by s ^ [0, /(S)] — )■ x{s) a length parameterization of Ti, there are Lipschitz 
parameterizations s € [0, /(S)] for S" such that x"(s) tends uniformly to x{s), and 

the left and right unit tangent vectors of x"(s) tend uniformly to the left and right tangent 
vectors of x{s). 

Remark 2.1. When /(S) = 0, these conditions are reduced to the uniform convergence of 
x"(s) toward x{s). This convergence can be defined by a family of neighborhoods around 
each element of Cp' (5"^), which therefore is a topological space. 

Definition 2.4. A crown S : (A,/u) — )• Cp'^{S^) is a continuous and monotone map E from 
the interval {X,n) into Cp'^{S^) endowed with the partial order ^. If the map is defined on 
the closed/open interval we talk about closed/open crown. 

Definition 2.5. When the crown is closed and increasing, S(A) = H'^ is called the interior 
curve of the crown and = S'^ its exterior curve of the crown. If the crown is 

decreasing, these names exchange. The range of the open/closed crown is denoted by S(]A,/^[), 
S([A,/i]), respectively. The crown itself as an ordered family of curves will also be denoted by 

Remark 2.2. A (closed) crown is a homeomorphism from a closed interval [A,/i] to S[A,^], 
since the map S is continuous, and the interval [A,^] is compact. 

Definition 2.6. Let A C ^ be a bounded, connected set, whose boundary consists of a finite 
number of disjoint Jordan curves. We call exterior curve of A the unique Jordan curve 
whose interior contains A and interior curve (s) the other Jordan curves S*, i € /. 

2.2. The class of Very Simple Functions. A digital image is usually known by its samples 
{uii, j)}o<,i<M o<j<N ^ rectangular grid of J] = [0, M] x [0, A^]. We assume that the 
underlying image < u{x) < 1 whose samples are the u{i,j) is a Lipschitz function defined 
on O, the continuous image domain. We shall always assume that u{x) = on the boundary 
of the domain dQ and that u{x) > in the interior of ft. 

The bilinear interpolation in il. is the simplest continuous interpolation from the discrete 
samples u{i,j). This interpolate, still denoted u, is defined as the unique function coinciding 
with the digital image u on the samples which is bilinear in each dual pixel (the square formed 
by the centers of four adjacent pixels). This means that u has the form 

u{xi,X2) = axi + (3x2 + 7x1X2 + S 

on each square with vertices {i,j), (i+l,j), («,j + l), (z+l,j + l). This bilinear interpolation is 
therefore positive on the interior of the domain and zero on dil.. The set of bilinear interpolates 
of digital images on il. will be denoted by BC{il.). The next result is a sane consistency property 
of the bilinear interpolation. 
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(a) (b) (c) 

Figure 1. Topological structure of the family of level lines, for a bilinear interpo- 
lated image. For every level < A < 1, except for a finite set, the iso- level sets 
{u = A} are made of a finite set of piecewise Jordan curves. At critical levels 
Afc, they can have a rather complicated form: they can have T-junctions or reduce 
to segments (a), contain square pixels (thus being flat areas) or have self-crossings 
at saddle points (b). However, by fattening the iso- level sets corresponding to the 
critical levels, the topology of level sets becomes very simple: only Jordan curves or 
flat sets (c). 

Proposition 2.1. If u defined on Q is Lipschitz and only known by its samples, its bilinear 
interpolate converges uniformly to u when the grid mesh tends to zero. 

The references |17j . |16j . [5] show that the bilinear interpolation brings a long list of nu- 
merically useful topological properties. 

Proposition 2.2. (Properties of level lines of bilinear interpolates) 

(LI) For every level < A < 1 except for a finite set of levels Ai, . . . , A„ called critical, 

the iso-level set {u = A} is the disjoint union of a finite set of piecewise-C^ Jordan 

curves, denoted by {T,^'^)i^j^ where I\ is a finite set of indices; 
(L2) The open set Vl\u'~^{{\i, ... ,\n}) has a finite number of connected components. Each 

connected component is the range of an open crown t}^''^^ where v € {Ai, . . . , A„}. 

As a consequence, is partitioned in open crowns and in the closed iso-level sets 

{n = Aj} corresponding to the critical levels. 

Sketch of proof. A dual pixel will contain a critical level either if it is flat, or if it contains a 
saddle point. In both cases there is only one critical level in the pixel. Since there is a finite 
number of pixels, there is a finite number of critical levels. At any other level, the restriction 
of an iso-level set to a given pixel is either empty, or is a single piece of hyperbola. Since 
the bilinear interpolate is continuous, these pieces of hyperbolae concatenate at each pixel 
boundary to form one or several disjoint Jordan curves. □ 

The above structure is quite simple, but it does not describe the structure of level lines 
at the critical levels Afe , which can take a rather complicated form as illustrated in Figure [TJ 
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To define a curvature evolution for these sets would be cumbersome. We shall overcome this 
drawback, both numerically and analytically, by building a still simpler approximation. This 
approximation of the image is obtained by fattening all critical iso-level sets into open sets 
whose boundary is a finite set of Jordan level curves. Their curvature evolution will simply 
be defined by the curvature evolution of their boundary. Thus, we define a still simpler 
approximation for Lipschitz functions, deduced from the bilinear interpolation. 

Definition 2.7. We say that a Lipschitz function u on 0, is very simple and we denote by 
u € ViS(rj) if it satisfies properties (LI), (L2) and 

(L3) Each crown can be completed into a closed crown by adding its interior and 

exterior curves S'^ and Yj^ , which are limits in the Cp' (S^) topology of the level lines 
of the crown. 

Remark 2.3. The construction of the partition goes as follows. The domain O is the union of 
all its isolevel sets {u = A}. By the property (LI), these sets are disjoint unions of (smooth) 
Jordan curves, except for a finite number. Property (L2) further tell us that these Jordan 
curves are organized in a finite number of open crowns. At this level, is the union between a 
finite number of open crowns and the finite number of iso-level sets leftover (which in turn can 
be split in a finite number of connected components, since we lay on a compact). The third 
property (L3) tells us that the boundaries of these flat areas, which are not necessarily Jordan 
curves, can be shifted to the crowns, such that the crowns become closed. As a consequence, 
for very simple functions, is the disjoint union of the ranges of a finite number of compact 
crowns, and of a finite number of flat open connected components, each belonging to some 
critical level A^, 



where K is the finite set of all crown indexes, fJ,k,^k ^ {^ii " " " j-^n} are critical levels, and 
Jj is the finite set of indexes of open connected components of the iso-level set of u at the 
critical level Xj. 

Definition 2.8. We call flat regions of uq € V5(r2) the open sets F'^i'\ whose boundaries 
are unions of a finite number of piecewise-C^ Jordan level curves, and on which the function 
is constant: 



Lemma 2.1. Every Lipschitz function can be approximated uniformly by a sequence of very 
simple functions. 

Proof. For all u € Lip{0,) and e > we build a function v € V5(r2) such that 





u{x) = Xj, Vx G F^^'\ I G Jj. 



sup \u{x) — v{x)\ < e. 
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By Proposition 12. 11 we can assume that u G BC{Q). We can order the set of its critical levels 

Xj\, the function v{x) = /^(^(x)), 



Ai < A2 < ... < A„. Then define, for < e < | min^ |Aj+i 
where / is the 1-Lipschitz nondecreasing function 

' t, t<Xi-e 

Xk-ke Xk- e <t< Xk + e, 

t — ke Xk+ e <t < Afc+i — e 

t — ne, t>Xn + s 



(2.6) 



fe{t) 



k 
k 



1, n 



l,n- 1 



We have — t\< ne, thus max^ |ii — f | < ne which proves the convergence claim. 

It remains to check that the function v € VS{Q). The critical levels of v are inherited from 
u and consist of {A^ — A;e}fc=i By construction, all the level lines of v at noncritical levels 
are level lines of u at noncritical levels and hence Jordan level curves, grouped in crowns. 
Consider now the flat regions of which are the open connected components of the sets 

{Ajt — e < u < Ajt + e}. A; = 1, n. 

Since e < ^minfc(Afc+i — A^), the boundary of each one of these flat regions is contained in 
the union of the level sets {u = A^ — e} and u = {Xt + e} which are a finite number of Jordan 
curves. Each one of these Jordan curves belongs to a crown Y^'^'^^k'i of u which is truncated 
into a crown Y}^k~ke,\y-k'e[ q£ ^_ □ 



3. Level Lines Shortening operator 



3.1. LLS semigroup operator. Given a very simple function uq € VS{^), its Level Lines 
Shortening evolution consists in evolving independently and simultaneously by Curve Short- 
ening each of its level lines, denoted by Sq'*, and eventually reconstructing, for each time 
t > 0, a new function u{-,t) whose level lines are the evolutions S^'*, where the subscript t 
denotes the time t. This definition, will have to be proven consistent. Our goal is therefore 
to prove the commutative diagram: 



MCM/LLS 

\ 

U{;t) 



level lines extraction 



reconstruction 



CS 



{^t''}x,- 



To this end, we shall use several fine properties of curve shortening evolution 
is given in terms of a nonlinear geometric partial differential equation 

dx 



13], which 



(3.7) 



dt 



{s,t) = k{s,t) 



where x{s,t) is a family of smooth Jordan curves parameterized for each t by a length pa- 
rameter. The vector k is the acceleration which is normal to the curve, points towards its 
concavity, and whose norm is the inverse of the radius of the osculating circle. 



Theorem 1. The curve shortening evolution has the following properties: 
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(PI) For any Jordan C^'^{S^) curve Sq, there exists a collapsing time T(So) > such that 
the Cauchy problem {3. 7p has a unique solution Sj in 

cO'i (51 X [0, r(So))) n [s^ x (o, r(So))) 

which still is a Jordan curve. We set T,t = 9 for t > T(So). 

(P2) The map Sq i— )• T,t is continuous for the Cp'^{S^) topology of Definition [KM For 
t = T(Y,q) the curve collapses to a point x{T,q). 

(P3) Before collapsing at time T{T,q), the curve T,t - rescaled at constant area equal to vr 
- converges in the Cp' (S^) topology to the unit circle centered at the collapsing point 
x(So). 

(P4:) Inclusion Principle: 

• i/ ^ S^, then J^j ^ T.} for all t > 0. 

• i/S^ ^ S2, then T^j ^ Y^f for all t > 0. 

(P5) The min-distance of any two disjoint curves increases with time until one of the curve 
collapses 

dist{T,l,T,'i) < dist{El,T,f),ys < t. 
(P6) Convex curves remain convex and shrink in time: Tit ^ Tq. 



Theorem 2. Properties (PI) — (-P6) hold for affine shortening: 
(3.8) ^(s,t) = |k|-2/3k(s,t) 

except for {P2>) which is replaced by 

(P3)' Before collapsing at time T(So), the curve Sf converges in the Cp'^{S^) topology to 
an ellipse centered at the collapsing point x(Eo). 

Let uq € VS{Q) be a very simple function and {Sq''}^ jg/^ its level lines. Denote by S^'* 
the evolution of Sq'* at time t 

Our first purpose is to show that the family of smooth Jordan curves {S^'*};^ jg/^ is actually 
the set of level lines of a very simple image u{-,t). 

Definition 3.1. Let Sq'"''^' = (Sg)^^^^^^ be a closed crown. We call level lines shortening 
of the crown Sq^''^^ the family of curves 

LLS{t) := (S,^)Ae[C,M]- 

where are the curve (affine) shortening evolutions ofT^ for all A € [C)A*]- 

To fix ideas we refer in the following to increasing crowns. Analogous results hold for 
decreasing crowns. 
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Proposition 3.1. Consider a closed increasing crown S'^'^l. Then the collapsing time T(A) = 
T{T,^) of the curves of the crown is a continuous increasing function of X & [C^lA- level 
lines shortening at time t < T(S'^) transforms S^^'^] into a closed crown 

LLS{t) (^C.A^l) =5][--(C,T-W),,]_ 



Proof. Since the composition of two continuous maps is continuous, and the composition 
of two strictly monotone maps is strictly monotone, this is an immediate consequence of 
Theorem [1] and of the definition of crowns. By property (-P1) level lines shortening preserves 
space-time continuity, whereas by (-P4) preserves strict monotonicity. 

□ 

In short, a crown remains a crown by level lines shortening, and is made of all curves of 
the initial crown which have not collapsed yet. It is convenient to also define the evolution of 
a flat region. 

Definition 3.2. Let Fq-''\ I £ Jj be a flat region of uq € VS{Q) at level Xj, of exterior curve 
T,Q and interior curves Yl^, m £ M. We call the level lines shortening of the fiat region 

Fq^' , I ^ Jj the set defined by 

LLS{t){F^^'') = F,'-' := /nt(Sf) n (^f]Ext{J:T)^ , Vt < T(S^) 

where S| and SJ", m (z M are the curve shortening evolutions o/Sq and S™, m £ M. 

The initial flat region remains a region whose boundary is made of all Jordan curves of 
its initial boundary which have not collapsed. By the inclusion principle the last curve to 
disappear is the external boundary. When it collapses, the region disappears. 

Theorem 3 (Definition of LLS for very simple functions). Let uq be a very simple Lipschitz 
function, with critical levels {Xk}k=i,- - ,n, Jordan curves Sq'* indexed by their level X and 
i £ L\, and of flat regions F^^'^ at critical levels Xj, indexed by I E Jj. The level lines 
shortening evolution of the function uq is the function LLS{t){uQ) = u(-,t) defined by 



A, if x £ S^'* 
Xj, ifxeF^"^ 



(3.9) u{x,t) 

'^0, ifxen\nt 

where fit is the domain surrounded by the curve shortening evolution {d^})t of the domain 
boundary dQo (which is the only zero-level curve of uq). Then this definition is complete, 
consistent, the evolved function is a very simple function u{-,t) S VS{Q) whose Lipschitz 
constant is smaller than or equal to the initial one. 

Proof. The initial domain is partitioned in crowns and flat regions whose boundaries are 
either interior or exterior curves of crowns, or dQ. By the min-distance property {P5) in 
Theorem [H when time increases the evolved level curves of uq fall apart from each other and 
so do the boundary curves of the flat regions. Thus, by Proposition [1] and Definition 13.21 
the crowns never meet and the flat regions are at all times the connected components of the 
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(a) Image associated to a very simple function and its corresponding level lines. 




(b) The LLS operator maps very simple functions onto very simple functions. 



Figure 2. The Level Lines Shortening (LLS) numerical chain: from the original 
image (top- left), the family of all its simple level lines is extracted at quantized non- 
degenerate levels (top-right). Simultaneously and independently, each level line is 
evolved (bottom-right) and the evolved image having these level lines is reconstructed 
(bottom- left). Both the images and their families of level lines satisfy the topological 
properties (LI), {L2), {L3). Performing the LLS evolution, digitization artifacts due 
to noise, compression and under-sampling are attenuated. 



complement in of the union of crowns. In other terms the evolved crowns and evolved flat 
regions form a partition of Qt given by 



(3.10) 




u 



u 



je{i,-,n}, 
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On the other hand, the boundary of $7 is convex and remains convex by curve shortening 
(property (-P6) in Theorem [1]). Hence if points x initially belonging to 17 have been crossed 
by the evolving boundary, we have u(x, t) = 0. This renders the definition complete and 
consistent. 

Let us show that u{-,t) is Lipschitz. A function is L-Lipschitz if and only if the min-distance 
between any two of its level sets with levels A and is larger than ^—j^ ■ More precisely, since 
the min distance between a flat region and another, or between a flat region and a level line, 
is always attained on the level lines which bound the flat regions, we have 

\C-f^\ 

' dist{Yl'^, Sq) 

Since all min-distances between level lines increase (property (-P5), Theorem [1]), we have 



Lip{uo) = max{- — - — ^ — —}. 



\C-f^\ ^ IC-^I 



and hence the Lipschitz constant of u{-,t) is smaller than or equal to the Lipschitz constant 
of Uq. □ 

Definition 3.3. We call Level Lines Shortening operator, shortly LLS, the above oper- 
ator acting on the class of very simple functions, 

LLS{t) : VS{n) ^ VS{Q) 

Uq u{-,t). 

Since the curve shortening itself is a semigroup, LLS{t) also is a semigroup, namely 

LLS{t + s)uo = LLS{t){LLS{s)uo). 
Corollary 3.1. The level lines affine shortening evolution of a very simple function 



(3.11) LLAS{t){uo){x) 



A, if X ^ S^'* 
A,-, ifxeF^^^'' 
0, ifxe^}\^}t 



where S^'* and F^^'^ are the affine shortening evolutions of the initial level lines, respectively 
flat areas, is well defined and maps VS{il.) onto itself, preserving the semigroup property. 

Proof. This comes immediately from Theorem [2l which ensures that all the topological prop- 
erties of curve shortening hold as well for affine shortening. □ 



3.2. Properties of Level Lines Shortening. 

Lemma 3.1. Let <j)Q be a radial increasing function centered at xq, i.e. (po{x) = <f{\x\) with 
if increasing. Then its level lines shortening evolution is given by 

<P{x, t) = ip (V|x|2 + 2(t-to)) 
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and its affine shortening evolution is given by 

Hx,t) = ^ (^(^\x\l + ^{t - to)) 

Lemma 3.2 (Local comparison). Let uq be a very simple function, (pQ a radial increasing 
function centered at xq, and denote by u{-,t) and (j){-,t) their LLS/LLAS evolutions. If in a 
local neighborhood M{xq) of xq the following holds 

uo{x) < (j)o{x),yx e M{xo) 

then there exists a short time to > such that 

u{xo,t) < ip{xo,t),yO <t <to. 

Proof. We argue for LLS, the case LLAS being analogue. The level set {i^o < ^} is the open 
disk with radius (j)Q^{X) and satisfies the inclusion 

{x G A/'(xo); (l)o{x) < \} C {x e Af{xo);uo{x) < A} 

Then every level line S^'*^'^ of surrounds no point belonging to a level set of the same level 
of uq. By the inclusion principle and the topological structure of the level lines for a very 
simple image, this property is preserved for all t < to where to is the vanishing time for the 
largest level line of (p in the neighborhood J\f{xQ). Indeed if the level set of u{x,t) is a flat 
part bounded by Jordan curves, by the min-distance property (-P5) these Jordan curves never 
cross in their evolution the circles corresponding to the level set for (j){x,t). If the level set 
of u{x, to) is a finite set of level lines, in the same way the evolved level lines of u{x, t) never 
cross the circular level line of (j){x,t). 

On the other hand xq belongs to all the level lines of </>(•, t), for all times t < to- Hence 
the value of u{-,t) at xq must necessarily be less than the minimum level of (p{-,t), which is 
attained exactly at xq. Consequently 

uixo,t) < (j){xo,t),yt < to. 

□ 

Proposition 3.2 (Space-Time Continuity). Let uo G V5(r2) be L-Lipschitz continuous and 
consider its level lines (affine) shortening evolution u{-,t) = LL{A)S{t){uo), for all t G (0, cxd). 
Then u G C°{n x [0,oo)). 

Proof. We want to find a Lipschitz type estimate for the function u and we argue separately 
in space and time: 

\u{x,t) - n(xo,to)| < \u{x,t) - u{xo,t)\ + \u{xo,t) - u{xo,to)\. 

By Theorem [3] the LLS evolution u{-,t) at any time t > of the initial function uq remains 
L — Lipschitz continuous and hence 

|ti(x, t) — u{xo, t)\ < L\x — xo\. 

The time-continuity follows from comparisons with shrinking cones. More precisely, the Lip- 
schitz continuity at time to tells us there exists </)(•, to) a L-Lipschitz radial upper barrier 
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u{x,tQ) < (j)(x,tQ), touching u at point xq, given by 

<j){x, to) = n(xo, to) + L\x - xo\. 
By Lemma |3.H its LLS evolution is 

(pix, t) = u{xo, to) + L-s/lx - xop + 2{t - to). 
It follows from the local comparison with radial functions given in Lemma 13.21 that 

u{xo, t) < (l){xo,t) = u{xq, to) + Ly/2{t - to). 

Consequently 

\u{x,t) - u{xo,to)\ < l(^\x - xo\ + \/2(t - to)) • 
This proves that u{x, t) is uniformly continuous in t and x. Similar results hold for LLAS. □ 

The structure of very simple functions implies that there is only a finite number of possible 
collapsing times, namely times > where some crown collapses to a point. At these 
times tl < t1 ■ ■ ■ < there is one (or several) collapse pairs (t^\ S^*^'**^). The next lemma 
gives a stability property for flat regions. 

Lemma 3.3 (Flatness). Let uo be a very simple function and u{-,t) its LL(A)S evolution at 
time t. Let xq be a point in a flat region of u{-,to) and suppose that it is not a collapsing 
point. Then there exists (5o > such that xo stays in aflat region ofu{-,t), for all |t — to| < ^0- 

Proof. Since xo belongs to a flat area of u{-,to) there exists a small ball B{xo,ro) centered at 
Xo meeting no other level line of u{-,to): 

(3.12) S,Yni?(xo,ro) = 0. 

The number of collapsing points being finite, we can also choose tq small enough, so that 
B{xo,ro) contains no collapsing point of the evolution of u. Let dB{xo,r{t)) be the circle 
centered at xo and evolving by Curve Shortening such that r{to) = rQ. 

1. Fix 5i = r^/A such that dB{xo,r{t)) has not collapsed at time t = to + ^i- Then it 
follows from (|3.12|) and from the inclusion principle that for all t < to + 6i no level line of 
n(-,t) meets the ball B{xo,r{t)) 

E^'ni?(xo,r(t)) = 0, 

where r{t) = y^rg - 2(t - to). 

2. We now prove that there exists a time §2 and a radius p > such that for all t G 
(to — '^2, to) no level curve of u{-,t) meets B{xo,p). Assume by contradiction that there exist 
tj, Xj and curves S^;' of LLS evolution of u such that 

(3.13) tj — )• to, Xj Xo with Xj G "^tj • 

Consider the corresponding initial curves Sq-' . Since the ball B{xo,ro) contains no collapsing 
points of u, the evolutions S^-' will have uniformly bounded curvature, and thus stay uniformly 
bounded in the Cp'^{S^) topology. It is then possible to extract a converging subsequence to 
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Sq. By the continuity property {P2) in Theorem [U it follows that T,^° contains xq, which 
contradicts our initial assumption. 

□ 



4. Equivalence with the curvature motions 

4.1. Mean Curvature Motion. We rewrite the geometric curve shortening in the form 
(CS) I = k., 

where k denotes the scalar curvature and z/ the unit normal to the curve, with i' continuous 
and the sign of k guaranteeing that ki' points towards the interior of the domain surrounded 
by the curve at convex points and towards the exterior at concave points. We are interested 
in its equivalence with 

{ut = \Du\curv{u), in x [0,oo) 
u{-,0)=uo, on ]R2. 



(MCM) 
where 



n 2 

\Du\curv(u) = \Du\div{-——) = } [dij - )ux,x,- 

\Du\ ^-^ \Uur ■' 

i,i=i 



We refer to a viscosity solution for the parabolic PDE, which is defined in terms of point- wise 
behavior with respect to a smooth test function. We use herein the definition presented by 
Morel and Guichard in [15], which was proven by Barles and Georgelin [2] to be equivalent 
with the viscosity solutions given by Evans and Spruck in [1^ and Chen, Giga and Goto in 

m- 

Definition 4.1. A function u G C(M^ x [0, oo))nL°°(M^ x [0, oo)) is a viscosity sub-solution 

of (MCM) iff for each (p € C°°(M^ x [0,oo)) such that u — (j) has a local maximum at (xo,to) 
we have 

(ptixQ^to) < \D(j)\curv{(j)){xo,to) if D(j){xo,to) 

</'t(a;o, io) < if D(j){xQ, to) = and D'^(j){xQ,tQ) = 0. 

A function u G C(M^ x [0, oo)) nL°°(M^ x [0,oo)) is a viscosity super- solution of (MCM) 
iff for each (p S C°°(M'^ x [0,oo)) such that u — (p has a local minimum at {xQ^to) we have 

(pt{xo,to) > \D(p\curv{(p){xo,to) if D(p{xQ,tQ) 

(ptixo,to) < if D(p{xQ, to) = and D'^(p{xo,to) = 0. 

Remark 4.1. The above definition can be further simplified [TJ]: replacing "local maximum 
(minimum)" with "strict local maximum (minimum)" one obtains an equivalent definition of 
viscosity solutions. Furthermore, it is enough to consider test functions of the form (p{x, t) = 
fix)+g{t). 
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(4.14) 



Theorem 4. Let uq G V5(il). Then the Level Lines Shortening evolution of the function uq, 

u{x, t) = LLS{t)uo{x),yx G ]R^ Vt G [0, oo) 
is a viscosity solution for (MCM) with initial data uq, 

Proof. It is sufficient to check that u{x,t) is a viscosity sub-sohition. Analogous assertions 
hold for viscosity super-solutions. Let (j) ^ C°°(M^ x [0, oo)) such that u — cj) has a strict local 
maximum at (xo,to)- Adding if necessary a constant, suppose that 

ti(xo,to) = 4'{xo,to) = A 
u{x,t) < (j){x,t),y{x,t) G V 

where y is a small neighborhood of (xq, to)- The proof is completed by the next three lemmas, 
where we distinguish two situations: either the point xq is inside a flat region of ^o) (lemma 
14. ip . or it belongs to some level line, singular or not of this function (lemma 14.31 in the case 
where D<j){tQ,XQ) = 0, and lemma [52] when the gradient is not zero). □ 

Lemma 4.1. Let xq be a point in a fiat area ofu{-,to). Then 

D(j){xo,to) = and (j)tixo,to) = 0. 

Proof. By Lemma 13.31 the function u is constant in a small neighborhood J\f{xQ,tQ). From 
the local maximum condition we deduce that the point (xo,to) is a local minimum for the 
test function <j) G C°°, which yields the conclusion. 

□ 

We consider now the case when xq belongs to a level line S^^'" of the function u{-,tQ). By the 
construction of LLS{t)uQ this level line is following the classical curve shortening. 

Lemma 4.2. Let xq belong to a level line S^j" of the function u{-,t(j). Let (j) be a smooth test 
function such that at the maximum point (xo,to) of u — (j) 

D<)){xo,to) ^0. 

Then <j) satisfies 

(j)t < \D(l)\div{-^^) at {xo,to). 

Proof. 1. The non-degeneracy condition Z)(/)(xo, to) 7^ and the regularity of the test function 
(j) imply by the implicit function theorem that the iso-level set 

= {xGf^;0(x,to) = A} 
is a smooth graph in a neighborhood of xq. A unit normal vector of S^j''^ at point xq is 

z^'^(xo,to) = TKI\(^o,to). 
\D(I)\ 

On the other hand, xq belongs to the smooth Jordan level line 5^^!)" '^^ ^- local 
maximum condition at point (xo,to) the two graphs Ti^f and S^j" are tangent at xo and 
do not intersect in a small neighborhood of the point. Therefore, the unit normal vectors of 
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x(t) 












A, It 



Figure 3. The level line '''' (in red) stays on the same side of the level line E 
(in blue). For t < tg, we consider the backwards locations x(t) of the point xq on the 
level lines E^'" as well as the intersections y{t) of the normal direction with the level 
lines E^'"^ of the test function. 



these curves coincide up to their sign. We set 

z^"(xo,to) = T^Tr(^o,io)- 
\D(t)\ 

Furthermore, for short times t € {to — 5, + 5) with 5 > small enough, the A level set of 
(pi-jt) denoted by 

sf'"^ = {x(£ nt;Hx,t) = A} 

remains a smooth graph in a neighborhood of xq. By the maximum condition (j4.14|) . S^'"^ 
stays on the same side of S^'" (see Figure HTT]) . 

To fix ideas, suppose that for points x close enough to S^j", u{x) > A in the interior 
domain bounded by the level line S^j". Then the A— level line of the test function T^'^^'^ lies 
locally outside the same domain. In addition, the normal vectors i^"(a;o,to) and ;/'^(xo,to) 
point inwards the interior domain bounded by the level line S^j". 

2. We consider now the backwards locations x{t) on the curves S^'" of the point xq G S^j", 
curves which represent the curve shortening evolutions of some level line Sq'**. Let v{x{t),t) 
be the unit inward vector at point xq of the level line € S^j''. The vector points in the 
direction of xq — x{t) and hence there exists d{x{t),t) such that 

(4.15) x(t) = xo- d{x{t), t)u{x{t),t) 

By the smoothness property {PI) of Theorem [H 

lim u(x(t),t) = z^"(xo,to)- 

Let y{t) be the intersection point of the outward normal direction —i'{x{t),t) (we take here 
into account that we have considered upper level sets, i.e. v'^ points inwards) with the level 
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curve S^'*^ (situated outside S^j"). Then there exists d{y{t),t) such that 

(4.16) y{t) = xo-d{y{t),t)u{x{t),t) 

and D(j){y{t),t) ^ 0, since S^''^ remains a local graph for short times. From the implicit 
function theorem we deduce that y{t) is uniquely defined and varies smoothly in time. Thus 
the following limit exists: 

hm = ^(to). 

Furthermore, since S^'"'^ stays on the same side of S^'" we have 

d{x{t),t) <d{y{t),t). 

Taking the inner product with u{x{t, )t) in equations (|4.15|) and (|4.16|) and dividing by t < to, 
the previous inequality implies that 

{—. T ,l^{xit),t)) < {— ,u{x{t),t)). 

t — to t — iQ 

Passing to the limits as t — >■ to we have 

(4.17) (^(to),^-(^o,to)) < (||(to),z."(xo,to)) 
and taking into account that (CS) gives 

(4.18) (— (to),z."(xo,to)) = P(xo,to) 
we get 

fe"<(^(to),z^"(xo,to)). 

But and are ordered by inclusion and meet at point xq. Thus, their curvatures at xq 
are ordered k'^{xQ,to) < /c"(xo, to)- On the other hand, from the regularity of the test function 
and the fact that Z)(/>(xo,to) 7^ the curvature at xo can be expressed as 

fc'^(xo,to) = -div(-^)(xo,to) 
\D(t)\ 

Consequently 

(4.19) _div(^)(xo,to) < (^(to),i^"(xo,to)). 



3. The sequence of points y{t) found before belongs to the A— level set of the test function, 
thus we have 

</.(y(t),t) = A, for tG (to -5, to]. 
Differentiating this identity with respect to t one gets for t = to 

dy 

{-Q^{to),D(j){xo, to)) + (^t(xo, to) = 0. 
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But at the touching point xq the normal vector to S^J*^ can be expressed as 

-^(xo,io) = z^"(2;o,to)- 
\D<P\ 

Consequently the equality above at point (xo,to) becomes 

dy 

(t>t{xo,to) = -\D(j)\{xo,to){—{xo,to),lj''{xo,to)) 
which by inequality ()4.19p implies 



Mxo,to) < |i:>(/>|div(-^)(a;o,to))- 
\D(j)\ 



□ 



Remark 4.2. If a point xq is a collapsing point at time to, then 

u{x, t) = A, Vt € [to, to + 6) 
Therefore the test function satisfies D<j){xo,to) = 0. 

Lemma 4.3. Let xq belong to a level line S^j" of the function n(-, tg)- Let (f> be a test function 
such that D(f){xo,to) = 0, D'^(j){xQ,to) = 0. Then (p satisfies 0j(xo,to) ^ 0. 

Proof. By Remark 14. II we can assume that the test function has the form <j){x, t) = f{x) + g{t). 
By assumption, for every {x,t) in a neighborhood of {tQ,xo) we have 

(4.20) n(x, t) - fix) - g{t) < u(xo, to) - f{xo) - g{to). 

1. Assume first that S^j" is not reduced to a point and say it is the curve shortening 
evolution of the original level line Sq'". Denote by S^'" the intermediate evolutions for 
< t < to- Consider, as before, for short times t € (to — (5, to] the points x(t) belonging to 
S^'" such that x(t) xo as t — > to. Since n(x(t),t) = n(xo,to) from inequality (|4.20p we get 

(4.21) g{to)-g{t)<f{x{t))-f{xo). 
However, since the level line evolves by curve shortening, we also have 

(4.22) x(t) = xo + (t - to)/c(xo, to)i^(xo, to) + o(t - to). 

Substituting this asymptotic expansion in (|4.2ip and recalling that Df{xQ,to) = yields 

9'{to){to-t) <o{to-t), 
which implies ^'(to) < 0. Since g'(to) = 0j(xo,to), this proves the announced statement. 

2. The only case not treated by the above argument is when xo is the collapsing point of 
some level line Accordingly, to is its collapsing time. In this case, the level line does not 
have a normal direction at (xo,to) and consequently equation (j4.22p is not valid anymore. 

Nevertheless, by Theorem [H property (-P3), for short times t S (to — (5, to) the points x(t) 
lie approximatively on a circle of radius R(t) respectively, with -R(to) = 0. Thus 



(4.23) |x(t) - xol c:^ V2(to-t). 
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Substituting this asymptotic expansion in (|4.2ip and taking into account that 

Df{xo,to) = and D^f{xo,to) = 

we obtain 

g'ito)ito-t)<o{to-t), 

which imphes again g'(to) ^0. □ 

The Level Lines Shortening can be extended by density to an operator acting on the class 
of Lipschitz functions. 

Corollary 4.1. For each uq € Lip{il.) consider a uniform approximation by very simple 
functions Uq S VS{Q) and define its LLS evolution by 

LLS{t)uo{x) = lim {LLS{t)u}i{x)) , 

n— )-oo 

Then the right hand side member is well defined and is a solution of (MCM) with initial data 
Uq. 

Proof. From the previous theorem we know that the LLS evolutions of the very simple func- 
tions Uq 

= LLS{t)u'^ 

are solutions of the mean curvature equation. By the comparison principle we know that 
max(ti"'(x,t) - u^{x,t)) < max(uQ(x) - u^{x)). 

Since the sequence of very simple functions is uniformly convergent, we deduce that for each 
t > 0, the family {LLS{t)uQ}n C VS{il) is a Cauchy sequence in the || • ||oo norm. By the 
stability properties of viscosity solutions, the limit function 

u{x,t) = lim u^{x,t) 

n— >oo 

is also a viscosity solution of [MCM) and hence satisfies 

Lip{u{-,t)) < Lip^uo). 

□ 



4.2. AfRne Curvature Motion. Similarly, one can connect the affine shortening 
(AS) | = .V3, 

where v{-,t) is the inner unit normal vector of the curve x{-,t) and k{-,t) the signed scalar 
curvature corresponding the choice of with the affine curvature motion 

{ut = \Du\(curv(u)Y^^ , in x [0,oo) 
u(.,0)=no, on]R2. 
In the definition of these nonlinear evolutions, for x € M, x^^^ stands for sgn(x)|x|^/^. 
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Theorem 5. Let uq € V5(il). Then the Level Lines Affine Shortening evolution of the 
function uq, 

u{x,t) = LLAS{t)uo{x)yx G M^Vi G [0, oo) 
is a viscosity solution for (ACM) with initial data uq. 

Proof. The proof of Theorem H] is purely geometric, thus the same arguments apply for level 
lines affine shortening. 

1. When D(j){xQ,tQ) ^ we need to estimate from above the right hand side of 

dy 

M^Oiio) = -\D(p\{xo,to){—{xo,to),i'''{xo,to)). 
The proof is literally the same up to inequality ()4.17p 

{^{to),U^Xo,to)) < {^^{to),u^{xo,to))- 

The only difference it makes with the previous proof is when (AS) comes into play. More 
precisely the evolution equation (|4.18|) at the maximum point {xQ,to) should be replaced by 
an affine shortening evolution 

(4.24) {^{to), u-ixo,to)) = {k^{xo,to)f' 
On the other hand k'^{xo,to) < A;"(xo,to) which implies 

{k'^{xo,to)f^ = sgn{k'f')\k^{xo,to)\'/'' < sgn{k^)\k^{xo,to)\'/' = {k^{xo,to)f' 
But for the test function (j), k'^{xo,to) = —curv{(t)){xQ,to). Hence 

-{curv{<t>))'/^ < {^^{to),u^{xo,to)). 
from where we deduce the desired viscosity inequality. 

2. For the case D(j){xQ, to) = and D'^4>{xq, to) = it is sufficient to replace the asymptotic 
expansions ()4.22p by 

(4.25) x{t) = xo + {t- to) (A:(xo, to)) ^/^i/(xo, to) + o{t - to). 
respectively ()4.23p by 



4 



(4.26) \x{t)-xo\ ~ ^l(to-t). 

This concludes the proof. 

□ 



5. Numerical Implementation and Applications 

The discrete Level Lines Shortening Algorithm performs accurate sub-pixel evolution by 
mean curvature motion. The complexity of the algorithm is directly proportional to the total 
variation of the image, since it acts simultaneously and independently on all of the level lines 
of the image (of course up to a quantization step) . The main goal of the implementation is to 
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obtain and move level lines with arbitrarily high sub-pixel precision, overcoming thus all the 
drawbacks of finite difference schemes based on pixel approximations. Moving simultaneously 
level lines extracted with high sample precision allows straight level lines with high gradient 
to stand still with LLS, whereas they are diffused by FDS, even in its stack variant. Further 
example and details about the numerical implementation are given in [9]. The following 
example illustrates the recovery of shapes freed from their aliasing, JPEG, and noise artifacts. 

5.1. Image restoration and visualization. Aliasing due to pixelization is common in 
scanned documents. LLAS can be used for a graphic quality improvement smoothing contours, 
see Fig. [H After smoothing, pixelized level lines become accurate curves with sub-pixel control 
points, whose curvature can be faithfully computed. Thus the whole chain can be viewed as a 
numerical preprocessing before further numerical analysis and feature extraction. But there 
is also a strong interest in the direct visualization of the level lines and of the microscopic 
curvature map of an image. 




Figure 4. Top: original cartoon image and its corresponding bilinear level lines. 
Bottom: LLAS evolution and smoothed level lines. 



5.2. Fingerprints restoration and discrimination. Minutiae such as cores, bifurcations 
and ridge endings characterize uniquely fingerprints. Their detection requires a careful 
smoothing, particularly to avoid a spurious diffusion mixing the ridges. The main objective 
of smoothing is to sieve the curvature extrema, which allow the fingerprint discrimination. 
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Figure 5. Original fingerprint and its level lines, Level Lines AfRne Shortening of 
the image and the corresponding level lines. 
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